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ABSTRACT 

We investigate chaos synchronization between Josephson junctions coupled uni-directionally with 
time-delay. We demonstrate the possibility of high quality synchronization with numerical simu¬ 
lations of such systems. The results are of certain importance for obtaining the high power system 
of Josephson junctions, which is promising for practical applications. 
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I.INTRODUCTION 

Chaos synchronization [1] is of immense fundamental importance in a variety of complex physi¬ 
cal, chemical, power, biological, economical and social systems with possible application areas in 
secure communications, optimization of nonlinear system performance, modeling brain activity 
and pattern recognition phenomena, avoiding power black-out, obtaining high power radiation 
sources, etc. 

Finite signal transmission times, switching speeds and memory effects makes time-delay sys¬ 
tems ubiquitous in nature, technology and society [1-2]. Therefore the study of synchronization 
phenomena in such systems is of considerable practical importance, see, e.g.[l-3] and references 
there-in. 

Josepshon junctions are essentially nonlinear dynamical systems [4], Study of chaos and its 
control in such systems [4-5] could be of huge practical importance from the application point 
of view. Chaos control in such systems is important for Josephson junction devices like voltage 
standards, detectors, SQUID (Superconducting Quantum Interference Device), etc. where chaotic 
instabilities are not desired. Chaotic Josephson junctions can be used for the short-distance secure 
communications, ranging purposes, etc [6]. In a chaos-based secure communications a message 
is masked in the broadband chaotic output of the transmitter and synchronization between the 
transmitter and receiver systems is used to recover a transmitted message [1-3]. As for ranging, 
as demonstrated in [7] the high bandwidth and rapid decorrelation of the laser pulse train could 
provide an unambiguous autocorrelation signal for measuring round-trip delay time. Compared 
with commercial ladar systems which use the time-of-flight measurement chaotic system use en¬ 
ables better accuracy. 
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It is well-known that Josephson junction is a source of radiation with frequencies np to tera¬ 
hertz region [8]. Terahertz waves, which lie in the part of the electromagnetic spectrum between 
infrared and microwaves, can penetrate materials that block visible light.They have a wide range 
of possible uses, including chemical analysis, security scanning, remote sensing chemical signa¬ 
tures of explosives in unopened envelopes, non-invasive detection of cancers in the human body, 
and high-bandwidth telecommunications [9-10]. In connection with the telecommunications ap¬ 
plication, it is noted that the earth’s atmosphere is a strong absorber of terahertz radiation in 
specific water vapor absorption bands, so the range of terahertz radiation is limited enough to 
affect its usefulness in long-distance communications. However, at distances of 10 meters the 
band may still allow many useful applications in imaging and construction of high bandwidth 
wireless networking systems, especially indoor systems. In addition potential uses exist in high- 
altitude telecommunications, above altitudes where water vapor causes signal absorption: aircraft 
to satellite, or satellite to satellite. 

However, the radiation from a single Josephson junction is very weak. Synchronization of the 
arrays of Josephson junctions is a natural way to increase the radiation power such sources. Upon 
achieving synchronization the radiation power will be proportional to the number of Josephson 
junctions squared. Recently it was established that certain highly anisotropic cuprate super¬ 
conductors naturally realise a stack of strongly clectromagnetically coupled intrinsic Josephson 
junctions [8]. From the fundamental point of view synchronization of huge number of complex 
systems could allow for higher efficiency of computational handling of such systems at least along 
the synchronization manifolds [1-2], 

In this paper we investigate chaos synchronization between three resistively-capacitively(RC)- 
shunted Josephson junctions coupled uni-directionally with time delays. We demonstrate the 
possibility of high quality synchronization with numerical simulations of such coupled systems. 
These results could be of certain importance for obtaining high power system of Josepson junc¬ 
tions, which is promising for practical applications. 


II.CHAOS SYNCHRONIZATION BETWEEN UNIDIRECTIONALLY COUPLED 
RC-SHUNTED JOSEPHSON JUNCTIONS 


Consider synchronization between the following RC-shunted Josephson junctions [5,11] (Fig.l and 
Fig.2) written in the dimensionless form: 
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Where 01,02,and 0 3 are the phase differences of the superconducting order parameter across 
the junctions 1,2 and 3, respectively;/? is called the damping parameter ((3R) 2 = h(2e/ c C')~ 1 , 
where J c , R and C are the junctions’ critical current, the junction resistance, and capacitance, 
respectively;/? is related to McCumber parameter /? c by /? 2 /? c = 1; /* is Planck’s constant divided by 
27r;e is the electronic charg e]id c i,dc 2 ,dc 3 is the driving the junctions direct current hoi ^ 2,03 cos(fii i2i3 t+ 
6*12,3) is the driving ac (or rf) current with amplitudes *01,02,03, frequencies fli,2,3 and phases 60 2i3 ;ti 
and r -2 are the coupling delay times between the junctions 1 — 2 and 2 — 3, respectively;coupling 
between the junctions 1 — 2 and 2 — 3 is due to the currents flowing through the coupling resistors 
R s 1 (between the junctions 1 and 2) and i?. s2 (between the junctions 2 and 3);a s i jS2 = -60201, 2 -60102 
are the coupling strengths between junctions 1-2 and junctions 2-3. We note that in Eqs.(l- 
3)direct current and ac current amplitudes are normalized with respect to the critical currents for 
the relative Josephson junctions; ac current frequencies fli j2 , 3 are normalized with respect to the 
Josephson junction plasma frequency u 2 x J2 J3 = 2e/ c i iC2iC3 (hC'i )2i3 )~ 1 , and dimensionless time is 
normalized to the inverse plasma frequency. 

It is noted that if the external drive current is purely dc then we have a second order au¬ 
tonomous system, Eq.(l) and therefore chaotic dynamic for Eq. (1) is ruled out. In order to 
make the dynamics of the Josephson junction (1) chaotic, one can add ac current to the external 
driving, which makes Eq.(l) a second order non-autonomous system. Then treating QR term as 
a new dynamical variable Eq.(l) (and consequently Eqs.(2-3)) can be rewritten as a system of 
third-order ordinary differential equations, which can behave chaotically. 
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We underline that the value of phases # 12)3 for the ac driving forces can be incorporated into the 
initial conditions for Eqs.(4-12). 

We note that Eq.(l) is also used for modeling of a driven nonlinear pendulum, charge density 
waves with a torque and sinusoidal driving held [12-13]. We also mention that there are some 
other ways to make the Josephson dynamics chaotic without resorting to the ac driving. One way 
is to use resistively-capaeitivcly-inductively-shunted (RCL-shunted) Josephson junction [5]. Addi¬ 
tionally one can make the Josephson junction dynamics chaotic via time- delays in the system [11], 


III.NUMERICAL SIMULATIONS 

Now we demonstrate that three unidirectionally time delay coupled Josephson junctions can be 
synchronized. We use DDE software in Matlab 7 for numerical simulations. We study chaos 
synchronization between Eqs.(4-12) using the correlation coefficient C [14] 

= < < x >)(y(t + A t)-<y>)> 

yJ< (x(t)— < x >) 2 >< (y(t + At)— < y >) 2 > 

where x and y are the outputs of the interacting systems; the brackets< . > represent the time 
average; At is a time shift between the systems outputs: in our case A t = 0. This coefficient 
indicates the quality of synchronization: (7=1 means perfect synchronization. 

First we simulate Eqs.(4-12) for the identical set of systems’ parameters:^ = 0 2 = 03 = 0.25, i dc 1 = 
i<ic2 — idc3 — 0.3, ioi = (02 = *03 = 0.7, Di = 0.2 = = 0.6, 6\ = 0 2 = #3 = 0,a s i = a S 2 = 

0.45, Ti = r 2 = 0.15.Throught the paper initial conditions are set to be different for the junctions. 
Fig.3(a) demonstrates dynamics of variable \\)\ ;hg.3(b)pictures the error dynamics 02 — 0 i;( 7 i 2 is 
the correlation coefficient between 02 and 0i. Fig.3(c) represents the dependence between variables 
03 and 02 with C 23 being the correlation coefficient between 02 and 03. In fig.3(d) we present the 
error dynamics between the end junctions 03 — 0i, while C 23 indicates the correlation coefficient 
between 0i and 03. These simulation results underline the high quality synchronization between 
three identical Josephson junctions with different initial conditions. 

Next we study the effect of the coupling delay time on the synchronization quality. It is worth 
noting that synchronization quality is very sensitive to the value of the coupling delay. As demon¬ 
strated by the numerical simulations with increasing coupling delays between the Josephson junc¬ 
tions correlation coefficients decays rapidly.In fig.4 we present the dependence of the correlation 
coefficients Ci 2 , C 23 and Cj 3 on the time delay. It is noted that we only increase the values of the 
coupling delay times, but keep them equal: r\ — r 2 . 
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Fig.5 shows the results of numerical simulations of chaos synchronization between 9 coupled 
Josephson Junctions with T\ = T 2 = • • • = r§ = 0.001. and the other parameters are as in 
hg.3: error dynamics of ipg — ipi'.Cig is the correlation coefficient between ipi and Jpg. Finally 
we briefly dwell on the role of noise in the synchronization. Noise can both create and deteriorate 
the synchronization quality between the systems. It all depends on the intensity of noise and 
parameter mismatches between the systems [15]. Detailed study of the influence of noise on the 
synchronization quality will be reported in the future. 


IV.CONCLUSIONS 

We have numerically simulated chaos synchronization between three uni-directionally time-delay 
coupled Josephson junctions. We have demonstrated the possibility of high quality synchronization 
between such systems. The results are important for obtaining high power system of Josephson 
junctions, which is promising for practical applications. 
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FIGURE CAPTIONS 


FIG.1.Schematic view of the three unidirectionally coupled Josephson Junctions -JJ\, JJ 2 and 
JJ3. Coupling is realised via the resistors R s 1 and R s2 . 

FIG.2.Schematic set-up of the RC-shunted Josephson Junction subject to the external driving 

I ext • 

FIG.3.Numerical simulation of Eqs.(4-12):(a)time series of Vh;(b)error 0 2 — 0i dynamics; Cj 2 is 
the correlation coefficient between the Josephson Junctions JR and JJ 2 ;(c)dependence between 
variables ^3 and 0 2 : C 23 is the correlation coefficient between 0 2 and 0 3 ;(d) error dynamics of 
03 — ipi'Ci 3 is the correlation coefficient between 0 3 and 0i. The parameters are:/3i = /3 2 = 0 3 = 
0-25, idd = idc 2 = idc3 — 0.3, ioi = *02 = *03 = 0.7, fR = fl 2 = H 3 = 0.6, 9\ = 0 2 = 0 3 = 0,a s i = 
a s2 = 0.45,7 i = t 2 = 0.2. Dimensionless units. 

FIG.4.Numerical simulation of Eqs.(4-12):Dependence of the correlation coefficients C\ 2 (solid 
line), C 23 (dashed line) and C \ 3 (dotted line) on the coupling delay time between the junctions 
t 1 = t 2 . Other parameters are as in figure 3. Dimensionless units. 

FIG.5.Numerical simulation of chaos synchronization between 9 coupled Josephson Junctions with 
r 1 — t 2 — ■■■ — t s — 0.001. and the other parameters are as in fig.3: error dynamics of 0 9 —0i:Cj 9 
is the correlation coefficient between 0! and 0g. Dimensionless units. 
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